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We give a brief overview of systems that show spiral patterns and spatiotemporally chaotic states. 
We concentrate on two physical systems: (1) the oxidation of CO on Pt(llO) and (2) ventricular 
fibrillation in hearts. The equations that have been suggested as simple models for these two different 
systems are closely related for they are both excitable media. We present these equations and give 
a short summary of the phenomena they yield. 



I. INTRODUCTION 

Spiral waves are found in a wide variety of physical, chemical, and biological systems. These 
include convection patterns in cylindrical cells pi], patterns in the Belousov-Zhabotinsky chemical 



reactions the oscillations of cyclic- AMP in aggregating Dictyostelium Discoidium amoebae 

the oxidation of carbon monoxide (CO) on surfaces of platinum single crystals calcium 
waves in the cell cytoplasm of Xenophus oocytes, and in ventricular fibrillation |Q which is 
perhaps the most dangerous form of cardiac arrhythmia. A variety of mathematical models have 
been used to describe these phenomena; we refer the reader to the review by Cross and Hohenberg 
for a recent discussion of many such models The models range from cellular automata Q to 
deterministic partial differential equations, such as the Complex Ginzburg Landau equation (for 
oscillatory chemical reactions near onset) and equations in the FitzHugh-Nagumo class (for CO 
oxidation and ventricular fibrillation). We concentrate on the Fitzlfugh-Nagumo types of equations 
here. Our description is not very technical since the talk on which this article is based was aimed 
at an audience comprising scientists from different disciplines. We begin with a brief overview of 
phenomena followed by an introduction to equations in the FitzHugh-Nagumo class. We then present 
some results for the specific models suggested for the oxidation of CO on Pt(ffO) and ventricular 
fibrillation including our recent studies | pT[ . 

Many chemical reactions display spontaneous spatial self-organisation if the reactors are not 
stirred and the reactants are fed at a constant rate. Among these reactions the oxidation of CO on 
Pt(flO) has been studied in great detail The experimental setup consists of a cell containing 
a Pt wafer kept at a constant temperature. The reactants are fed into the cell at a constant rate 
and the products removed. Spatial structures in the concentration fields of the reactants have been 
observed by using photoelectron emission microscopy (PEEM). For a picture of the growth of a 
spiral wave on such a Pt(llO) wafer see Ref. [^. 

A human heart consists of four chambers - the two upper atria and two lower ventricles. Heart 
activity is driven by periodic signals emitted by the cells of the sinus node situated in the upper 
right atrium. These cells initiate the heartbeat cycle by emitting a pulse of electrical activity. 
This activity travels as a wave across the atria and stimulates the nerves of the atrioventricular 
node (AVN). From the AVN the signal travels along the His bundle and into the ventricles via 
the Purkinje fibre network which extensively innervates the ventricular wall. The ventricles now 
contract thus completing the beat cycle. Cardiac arrhythmia is caused when the signal from the 
sinus node does not propagate in the normal manner. Especially dangerous is ventricular fibrillation, 
an irregular pulsation of the ventricles, which is a major cause of death in the industrialised world 
|l^ . It has been conjectured for many years that ventricular fibrillation is associated with the 
formation of spiral waves (or their three-dimensional analogues called scroll waves) on the walls of 
the ventricles. Experimental evidence for such spiral waves has been increasing with the advent of 
voltage-sensitive dyes and advanced imaging techniques Figure 2 shows such a spiral wave on 
the ventricle of a canine heart. Such spirals have also been seen to break up. 
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FIG. 1. A spiral wave on a canine heart imaged using voltage-sensitive dyes and CCD cameras, after Ref. [ 
obtained from and reproduced with the permission of W. Ditto. 

Both the oxidation of CO on Pt (110) and ventricular fibrillation have been the focus of many 
theoretical studies. These are motivated partially by a desire to understand the phenomena seen in 
experiments and partially to elucidate the general properties of systems displaying such spatiotem- 
porally chaotic, statistical steady states. The remainder of this paper is organised as follows: Section 
2 consists of an introduction to FitzHugh-Nagumo equations and the models we consider. Section 3 
describes some results for the model of CO oxidation on Pt(llO). Section 4 describes some results 
for the model of ventricular fibrillation. Section 5 ends with concluding remarks. 



II. MODELS 



FitzHugh-Nagumo-type equations belong to the general class of reaction-diffusion equations. 
These are of the form 

f)TT. 

^ = /.(i7,)+A,V^[/,; 

where the fi are nonlinear functions of the Uj's and Dij are elements of a matrix of diffusion 
constants. The FitzHugh-Nagumo-type system used by Hildebrand et. al. for the oxidation of 
CO on Pt(llO) is 

du 2 1 / , N / r ,\,\ 

— = V M u(u — l)(u — (v + b)/a), 

here the fields u and v are related to the CO coverage and the surface reconstruction, a, b, and e 
are control parameters. Physically, e is proportional to the ratio of the rate constant for change in 
surface structure to the rate constant for the oxidation of adsorbed CO on the surface. Further, 
f{u) = if u < i , f{u) = 1 - 6.75u{u - 1)^ if | < u < 1, and f{u) = 1 if m > 1. The form of 
f{u) leads to the production of v only above a certain threshold value of u. This is an effective 
model with u and v related to the reactant coverages and the surface reconstruction. We refer the 
reader to Ref. for further details and the specific form for f{u), which has been found to model 



the CO oxidation on Pt(llO) experiments well. We have studied Model (Q) with both periodic and 
Neumann (no-flux) boundary conditions. Here we restrict ourselves to Neumann B.C.'s for ease of 
comparision with model (^). 

The Panfilov-Hogeweg model for ventricular fibrillation that we consider is |^ 

^=e{e,g){ke-g)- (3) 

where e is the transmembrane potential and g is a recovery variable corresponding to the amount 
of openness of the ion channels [^. Also, /(e) = Cie if e < ei; /(e) = — C2e + a if ei < e < e2; 
/(e) = C3(e — 1) when e > e2, and e{e,g) — ei when e < e2; e(e,(?) = £2 when e > e2, and 
e{e,g) = £3 when e < ei and g < gi. Physically e(e,g)~^ determines the refractory time period of 
this model of ventricular muscle. Since the ventricles are electrically insulated from the atria, we 
impose Neumann boundary conditions. 



III. CO OXIDATION ON PT(llO) 

For the model of CO oxidation on Pt(llO) (^, Hildebrand et. al. find that there are many 
different statistical steady states. These include states with no waves (N), states with flat shrinking 
waves (F), states with rigidly rotating spirals (S), states with meandering spirals (M) and turbulent 
states with creation and annihilation of spirals (Tl and T2). These authors have found a stability 
diagram or "phase diagram" in the b — e plane for a = 0.84, and have calculated the spatial 
autocorrelation function of the density of spiral cores {n{x, t)n{x + s,t + r)). They suggest that the 
change in the spiral dynamics on passing from Tl to T2 is like a "liquid-gas" transition. We have 
simulated model (^) by discretising it on a square lattice and using a variable stepsize fourth-fifth 
order Runge-Kutta integrator for time stepping j^. In some cases we also use a fast integration 
scheme proposed by Barkley We have used Neumann boundary conditions to compare with 

model (^). We have investigated the local dynamics, i.e., the dynamics obtained by plotting u{x'n,tn) 
versus v{x^, tn),n = 1,2,... for fixed spatial location xli and t„ = nA where A is a sampling interval 
which we choose to be = 0.1 time units. Such local phase portraits show that, in the states S, M, 
and Tl, the system displays oscillatory behaviour; a representative phase portrait is shown in Fig. 
(0) for Tl. 
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more disordered in the regimes M, Tl and T2. 
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The trajectory in the local phase portrait loops around the fixed point («, = (0.66, 0.484) of 
equation (p|) without the term. We define the phase ';/>(x,t) = tan~^ ((w(x, t) — ii*)/(it(x, t) —it*)) 
and note that it winds by 27r around the cores of spiral defects Thus it can be used to obtain 
the defect density p [ prf . 




FIG. 3. A plot of u{x) as a function of space in model (H). The colour at a point is proportional to the value of u at that 
point, obtained after evolving a typical initial condition for 5000 time units . Note the similarity of this figure with Fig 
For this figure e = 0.073 which corresponds to the Tl state. 




point is proportional to the value of the phase at that point. 



IV. VENTRICULAR FIBRILLATION 



It is believed that ventricular fibrillation occurs because of the spontaneous breakdown of contrac- 
tion waves in the heart muscle leading to a state with many rotating spirals Physiological 
abnormalities do not seem to play an important role in this arrhythmia in a sizeable fraction of 
cases. Thus it is interesting to study models in which spiral breakup appears spontaneously without 
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introducing spatial inhomogenieties in the diffusion constant, for example. For the model (g), Pan- 
filov et. al. M find that, for the parameter values ei — 0.0026, 62 = 0.837, Ci — 20, C2 = 3, C3 = 
15, a — 0.06, fc = 3,ei^^ — 75, £2^^ ~ ~ 1.8 and 0.5 < £3"^ < 10, an initial condition in 

the shape of a broken wave gives rise to states containing many spirals. They have also calculated 
an electrocardiogram (ECG) numerically for this model in the turbulent state. In this section we 
elucidate the spatiotemporal nature of the chaotic behaviour of this system. In Ref. |^ it has been 
noted that the breakdown of spirals in models ^ and 1^ is qualitatively similar. We elaborate on 
this similarity below. The state with broken spirals in model (^) is similar to the state Tl in model 
(H) as can be seen by comparing Figs. (Bl) and (H). 




condition for 1000 time units. 

We have simulated model (^) using finite-differencing in space with a forward euler algorithm for 
time stepping. Most of our simulations are done with a system size of 256 x 256 with a spatial grid 
step of da; = 0.5 and a time step of dt — 0.022. Physically, one time unit is equal to 5 milliseconds, 
while one spatial unit is equal to 1 millimeter. We use the parameter values given above. We have 
calculated a local phase portrait for this model, i.e. a plot of e{x„, tn) vs g{xn,t„) for a fixed spatial 
location Xn- This is shown in Fig. Mj. Note that this local phase portrait is similar to that in the 
Tl state of the CO oxidation mode02|). 
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FIG. 6. The upper figure shows a local phase portrait of the model 
potential e. The system size was 256 x 256 with a grid spacing of 0.5. 
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f) ; the lower shows a time series of the transmembrane 



We find that for a system of size 50 x 50 the lifetime rt of the turbulent state is ~ 300 time 
units. After this time the entire field configuration goes to zero for various initial conditions. For 
a system of size 100 x 100 Tt ~ 950 time units. We conjecture that rt grows with system size; we 
are investigating this scaling of rj. The typical size of a human heart corresponds to a system size 
of ~ 600 X 600 so the size dependence of rt is of obvious interest. We note in passing that this 
qualitative trend of rt increasing with system size is in accord with the well-known phenomenon 
that only the larger mammals have heart attacks. 

To characterise chaos in model (^), we have calculated the maximum Lyapunov exponent \m.ax 
for a system of size 128 x 128. The maximum Lyapunov exponent measures the divergence of nearby 
trajectories in the system. To calculate it we use standard methods . In the transient turbulent 
state of model (^), \max saturates to a value of ~ 0.2. This is shown in Fig (^). Note that when 
this transient decays, the maximum Lyapunov Exponent decreases to negative values. 
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FIG. 7. The approximation to the maximum Lyapunov exponent Xmax at time t versus t. Note that Xmax approaches a 
positive constant (~ 0.2) and then decays at large times to negative values. This occurs because there is a long-lived chaotic 
transient which finally decays to a state that is constant in space with e{x, t) = and g(x, t) — everywhere. The lifetime of 
this chaotic transient increases with the size of the system (see text). 



V. CONCLUSIONS 

Spiral patterns are seen in a wide variety of natural systems Q. Here we have concentrated on 
two such systems which are also excitable media. The chaotic states of these systems display the 
breakup of large spirals into smaller ones in the states Tl and T2. The qualitative similarity of such 
spiral breakup in models ^ and has been mentioned briefly earlier Q. Here we have elucidated 
this similarity by comparing the spiral patterns (Figs. ^ and 0) in the state Tl of model (^) and 
the chaotic state of model (^; local phase portraits (Figs, g and ^ are also similar. We have 
characterised the chaos in model (|^) by computing the maximum Lyapunov exponent Xmax (Fig. 
^) as we have done for model (^ elsewhere 

Our study shows that there is considerable similarity between the behaviour of model ^ in the 
state Tl and model (^ for the parameter values given above. It should be possible to exploit this 
similarity to gain new insights about these two excitable media (e.g., by doing experiments on the 
oxidation of CO on Pt(llO) which might be hard to do on a human heart). Of course some caution 
must be exercised while doing this for the two models are different in some respects: the chaotic 
state of model (|^ is a transient for the system sizes we have studied, whereas the analogous state 
for model (^) seems to be a statistical steady state; furthermore, if the terms are dropped in 
both models, then model (H) has an extra unstable fixed point [at (u, v) = (u*, u*)], i.e., the reaction 
kinetics are different in the two models. We hope our work will stimulate experimental studies that 
will try to elucidate the similarities and differences between spiral turbulence in the oxidation of 
CO on Pt(llO) and ventricular fibrillation. Since ventricular fibrillation in cardiac muscle is fatal 
|l^,^,^, any insight gained about it by studying similar systems is of great interest. 
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